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Stereocilia are actin protrusions with 
remarkably well-defined lengths and 
organization. A flurry of recent papers 
has reported multiple myosin motor pro- 
teins involved in regulating stereocilia 
structures by transporting actin-regu- 
latory cargo to the tips of stereocilia. 1 13 
In our recent paper, we show that two 
paralogous class 3 myosins — Myo3a and 
Myo3b — both transport the actin-regu- 
latory protein Espin 1 (Espl) to stereo- 
cilia and filopodia tips in a remarkably 
similar, albeit non-identical fashion. 1 
Here we present experimental and com- 
putational data that suggests that subtle 
differences between these two proteins' 
biophysical and biochemical properties 
can help us understand how these myosin 
species target and regulate the lengths of 
actin protrusions. 

MY03A has been implicated in late onset 
hearing loss in humans and mice. 14 Recent 
studies showed that myosin III A (Myo3a) 
targets the deafness-associated actin- 
binding protein Espl 15 to the tips of actin 
protrusions. 316 These two proteins syner- 
gistically interact such that Myo3 motor 
activity helps concentrate Espl at actin 
filament barbed-ends, while Espl WH2 
activity enhances actin polymerization, 
thereby increasing steady-state filament 
(and therefore actin protrusion) lengths. 3 
The late-onset hearing loss phenotype 
associated with mutations in MY03A led 
us to hypothesize that the MY03A paralog, 
MY03B, could be partially compensating 



for MY03A. One key difference between 
Myo3a and Myo3b is that Myo3b lacks an 
actin-binding site in its tail that is essen- 
tial for its ability to target the tips of actin 
protrusions. 17 Interestingly, Myo3b and 
Myo3a constructs lacking the actin-bind- 
ing tail domain (THDII) can still target 
the tips of actin protrusions in the pres- 
ence of Espl, and this rescue of tip-local- 
ization is dependent on the C-terminus 
actin-binding site in Espl. 1 

The observation that unlike Myo3a, 
Myo3b can only localize to the tips of ste- 
reocilia when associated with Espl, 1 led us 
to further investigate how these proteins 
behave when all three are coexpressed in a 
heterologous system. We found that when 
Myo3a, Myo3b, and Espl are coexpressed 
in COS7 cells, Myo3a invariably accumu- 
lated at the tips of filopodia in a steep tip- 
to-base gradient, while Myo3b localized to 
filopodia tips proximally to Myo3a with a 
much shallower tip-to-base gradient along 
the shafts of filopodia. We also found that 
inversely proportional amounts of the two 
myosins were found in the overlapping 
regions near the tips of filopodia, strongly 
suggestive of competition between the 
two proteins for localization to the same 
compartment (Fig. 1). 

In order to test whether the differen- 
tial myosin distributions were the result 
of differences in the tail vs. the motor 
actin-binding domains of the two iso- 
forms, we also performed cotransfections 
with Myo3aATHDII, Myo3b, and Espl, 
as well as Myo3aATHDII, Myo3a, and 
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Figure 1. Filopodia from COS7 cells co-expressing GFP-Myo3a and cherry-Myo3b consistently dis- 
play GFP-Myo3a accumulated at their extreme tips, while cherry-Myo3b consistently trails behind 
Myo3a with a relatively longer tip-to-base decay length (left column). GFP-Myo3a with and with- 
out the 3THDII actin-binding site (GFP-Myo3aATHDII) consistently accumulates at filopodia tips 
ahead of GFP-Myo3b (left and middle columns, respectively), while GFP-Myo3a does not exclude 
GFP-Myo3aATHDII (right column). The black and gray boxes are the motor domains of Myo3a and 
Myo3b, respectively. The red crescents represent the IQ domains. 



D m f 



Fixedrrvl < ^ ► < 'I' Noflux 

actin filament 



base z=0 n tip z=h 



Figure 2. Schematic description of the calculated model of myosin distributions. We consider 
freely diffusing (yellow), processively walking (green) and stalled (red) myosins, that convert into 
each other (red arrows) with defined reaction coefficients. At the protrusion tip we illustrate the 
two types of boundary conditions that we used; walking motors detach and diffuse freely from 
the tip, or become stalled and are carried away by the actin treadmilling. 



Espl (Fig. 1). We found that, like Myo3a, 
Myo3aATHDII excluded Myo3b from 
filopodia tips. In contrast, we found that 
Myo3aATHDII was not excluded from 
filopodia tips when transfected with 
Myo3a. Finally, we also observed that 
Myo3aATHDII appeared to have a longer 
decay length than Myo3a in its tip-to-base 
gradient distribution, somewhat similar to 
Myo3b (Fig. 1). 

To gain further insight into which fac- 
tors determine the distributions for the 
different Myo3:Espl complexes, we gen- 
erated a computational model based on 
continuum equations to describe the dis- 
tribution of myosin motors along a one- 
dimensional protrusion [4]. The motors 
are described to be in one of three states: 
freely diffusing (m^, moving ("walking") 
processively toward the tip (wj at velocity 
v , or bound to the actin but stalled and 



therefore moving toward the protrusion 
base (m ) at the treadmilling velocity of 
the actin (v ) . This model is schematically 
shown in Figure 2. 

The dynamical equations are given by: 

^ = DV 2 m f + k off m w + k offib m s - 

)m / (l - (m w + m s )) 

*jjjr = -Om ~ v a )Vm w + k offiSW (l + 
K(m w + m s )P)m w + k onsw m s + k on m f (l - 
(m w + m s )) - k off m w 

7^ = v a Vm s + k offtbw (l + K(m w + m s f)m w - 

k 0 n,sw^s ~ k on m f (l - (m w + m s )) 

where D is the diffusion coefficient of 
the free myosin, and K is the coefficient 



describing the process whereby walking 
myosins become stalled due to interaction 
with fellow actin-bound myosins (stalled 
or walking). This interaction term could 
be due to "traffic-jam" effects as well as 
biochemical events, such as Myo3a out- 
competing Myo3b in binding to Espl. The 
power (3 in this interaction term describes 
the cooperative nature of these processes. 
We chose the parameters in Table 1 below, 
while all the other parameters remained 
unchanged and listed below. We also take 
into account in Equation 1 the excluded- 
volume effects between actin-bound 
motors. 

The boundary conditions we chose are 
of a constant concentration of myosins at 
the base of the protrusion, where there is 
contact with the large reservoir of the cell 
cytoplasm, and of zero-flux at the closed 
tip. Walking myosins that reach the tip 
are assumed either to fall off the actin 
filaments and become freely diffusing, 
or to become "stalled," whereupon they 
remain bound to the actin core which is 
slowly moving back toward the base due 
to retrograde flow and actin treadmilling. 
We tested these two different choices, and 
found that the best overall agreement for 
all the different experiments was captured 
by the boundary condition of myosins 
detaching at the tip. 

The situation in the case of two dif- 
ferent Myo3 species shown in Figure 1 is 
more complicated, since we now have to 
consider two sets of equations of the form 
of Equation 1. We expect that the actin- 
binding affinity of Myo3b is lower than 
that of Myo3a due to the lack of an actin- 
binding tail in Myo3b. This should in 
theory decrease the overall actin-binding 
affinity of Myo3b. Furthermore, we also 
know from in vitro experiments that the 
motor domain of Myo3b has a lower affin- 
ity for actin as well. Given the ^50% reduc- 
tion in ATPase and ADP-release rates for 
Myo3b along with our kymograph data 
showing ^2x faster tip-directed velocities 
for Myo3a vs. Myo3b, 1 we defined the 
Myo3b tip-directed velocity to be about 
half that of Myo3a. We first checked if 
these differences alone are enough to 
explain the distributions shown in Figure 
1, and we found that we could not get any 
agreement in terms of the mutually exclu- 
sive localization of Myo3a vs. Myo3b. 
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From this we concluded that there 
must be some form of direct interac- 
tion (beyond the excluded volume effect) 
between the myosin species. The most 
obvious form of interaction to consider 
was that walking Myo3b could stall much 
more readily than Myo3a when the two 
proteins collide, perhaps by having its 
Espl cargo stolen by neighboring Myo3a 
molecules. We included this effect (Xterm 
in Equation 1), so that this term becomes 
(for the dynamics of the walking Myo3b): 

~ K 'OffsW 

(1 + K(m wa + Tn sa ) a )m w b 

where the second subscript denotes the 
type of Myo3. With this description we 
were able to reproduce the main qualita- 
tive characteristics of the distributions 
shown in Figure 1. Interestingly, since 
Myo3aATHDII is not excluded from filo- 
podia tips in the presence of Myo3a, our 
results predict that the Myo3a THDI has 
a higher affinity for Espl than the Myo3b 
THDI. Future experiments will be able to 
determine whether this prediction is cor- 
rect, which we believe to be likely, espe- 
cially considering that the Myo3a THDI 
sequence is significantly different from the 
Myo3b THDI sequence. 18 This demon- 
strates the utility of this theoretical mod- 
eling approach in providing insight that 
can guide future experiments. 

It should be noted that these distribu- 
tions are not precise fits. It's important to 
consider that all of the interactions between 
the motors have been lumped into the sin- 
gle TTterm in our analysis for simplicity's 
sake. Furthermore, even within the sim- 
plified assumptions made in this study, we 
have not exhausted all of the regimes of 
solutions possible within the confines of 
the current model. Nonetheless, we found 
that changing these key parameters for 
each myosin construct allowed us to qual- 
itatively reproduce the main features we 
observed in these distributions — namely, 
the longer tip-to-base decay length for 
Myo3b and Myo3aATHDII (controlled 
by reducing their k and by increasing 
their tendency to stall or fall off), and the 
exclusion of Myo3b from the filopodia 



Table 1. Parameters used for simulations 





v m (|jim/sec) 


k on (sec 1 concentration 1 ) 


K (concentration - 1 ) 




Myo3aAK 


0.075 


1 


0 


0 


Myo3bAK 


0.05 


0.5 


5000 


1 


Myo3aATHDII 


0.075 


0.5 


0 


1 



k = k „ = k „sw = 0.1 sec 1 , k sw = Isec^concentration 1 , k s = O.lsec^concentratiorr 1 , and the COn- 
orr offs off 'on 'on ' 

centration of free myosins in the bulk (protrusion base) was taken to be: mf = 0.0001 in arbitrary 
units of concentration (the final plots in Figure 1 give relative concentration with respect to the 
peak value). The treadmilling velocity of actin was taken to be: v a = 0.01 (|jim/sec). We convoluted 
all of our calculated distributions with a 50 nm wide Gaussian function to account for our finite 
optical resolution. 



tips in the presence of Myo3a (which is 
controlled by making Myo3a able to out- 
compete Myo3b for Espl cargo binding, 
thereby rendering Myo3b unable to walk 
toward the tip). Future experimental and 
computational studies will have to be done 
to attempt to dissect the details of the con- 
tributions of different interaction terms. 

It is intriguing that Myo3a elongates 
Espl filopodia about twice as much as 
Myo3b. Our results strongly suggest that 
the different biophysical and biochemical 
properties of these two paralogous motor 
proteins are directly related to their pro- 
portionally dissimilar filopodia elonga- 
tion activities. Specifically, we believe that 
the filopodia elongation activity of class 3 
myosins is directly proportional to their 
tip-directed velocities. The direct propor- 
tionality between the tip-directed veloci- 
ties of Myo3a vs. Myo3b and the lengths 
of Myo3a vs. Myo3b filopodia is strongly 
suggestive, and is inherently logical: A 
faster motor gives rise to a larger flux of 
cargo proteins to the tip, which should 
result in greater polymerization activity 6 
and/or a greater flux of actin monomers, 19 
depending on the context or cargo of 
these motors. For example, since Espl has 
both crosslinking activity as well as actin- 
polymerization activity, 3 ' 20 a larger flux 
of such a protein near actin barbed-ends 
could enable a faster rate of actin polymer- 
ization and/or a more highly cross-linked 
actin bundle, which in turn can produce 
a longer protrusion. 19 Furthermore, our 
earlier studies demonstrated that Myo3a 
N-terminal kinase activity modulates 
Myo3a motor and elongation activity, 
since autophosphorylation of the motor 
domain by the Myo3a kinase domain 



lowers the motor domain's ATPase rate. 16 ' 21 
Taken together, all of our observations 
point toward a direct relationship between 
the efficiency of motor-dependent tip-tar- 
geting of actin-regulatory cargo and actin 
protrusion lengths. 

In summary, through our combina- 
tion of experimental and computational 
approaches, we were able to dissect the 
underlying molecular components of 
these myosin: cargo dynamics. This work 
demonstrates the intricate nature of the 
interactions between molecular motors 
within cellular protrusions, and the use- 
fulness of theoretical modeling in deci- 
phering them. Specifically, we can begin 
to ask and answer detailed questions about 
the contribution of different aspects of 
myosin motility, (e.g., actin-binding and 
motor activity), to not only their biophysi- 
cal behavior, but also their physiological 
function. Indeed, these results strongly 
indicate that Myo3b can behave quite 
similarly to Myo3a, but they also show 
that Myo3b "falls short" of Myo3a activ- 
ity in terms of elongation activity. This 
finding is one plausible explanation for 
the late-onset phenotype associated with 
mutations in MY03A. Overall, these 
studies lend stronger support for further 
exploring the potential role for Myo3b in 
compensating for Myo3a, and how these 
two motor proteins interact in a physi- 
ological context. 
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